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Abstract 

From the published results of experimental research of ^^Fe(n,27) reaction car- 
ried out in Budapest, the values of probable densities of cascade intermediate levels 
with 1/2, 3/2 spin and radiative strength functions of cascade El and Ml transi- 
tions in ^^Fe compound nucleus have been determined. These results correspond 
to analogous data for other nuclei studied by us and contradict to predictions of 
conventional models. 



Introduction 

The modern experiment techniques permit one to measure parameters of most nuclear 
reactions with error about 10% or less. Meanwhile, specific characteristics of nucleus which 
are important for theory and practice such as level density and radiative strength functions 
of cascade transitions are derived from experiment with less accuracy. The principal 
cause of this is insufficient volume of experimental information used in fitting process 
or incorrectness of its algorithms. The level density and strength functions of cascade 
gamma-transitions below neutron binding energy are typical examples of problems to 
elucidate experiment results with high reliability levels. This is certainly when: 

1. Setting up the best method to find unknown quantities from experimental spectra 
if they cannot be determined directly. 

2. Collecting the experimental data as much as possible to guarantee the singular- 
ity of the desired parameters and the ability to use these quantities for complementary 
independent revision. 

3. Evaluating all the sources of systematic error of achievable values. 
Actually, the level density p was determined now from experimental data on: 

(a) the nucleon evaporation spectra (see, for instance, [1]); 

(b) the intensities of two-step gamma- cascades [2] 

3;^^ Fa 'j;^ < T xi > mxi <Tif>mif 

between compound state (neutron resonance) and a group of low-lying levels of studied 
nucleus determined according to algorithm [3] for all possible intervals AE of energy Ei 
of cascade primary gamma-transitions; 

(c) the gamma spectra depopulating levels with the excitation energy Eex in (d,p) [4] 
and (^He,Q;) [5] nuclear reactions;. 



Within methods [2] and [5] , p is determined simultaneously with strength functions of 
cascade gamma-transitions: 



Here is the gamma-transition energy and Dx is the distance between initial levels A. 
Such form of k provides minimal dependence of this parameter on nuclear mass and allows 
direct comparison of these parameters in nucleus with different masses. It is assumed 
that the partial radiative width Txi of gamma transition between A and I levels is the 
random average value, and the width ratio for different multipolarity transitions can be 
determined experimentally. It is evident that the criterion of the necessary confidence for 
the determination of k and p does not cause the difference among the results of various 
experiments in principle. In other words, once there are errors, explanation for these is 
obligatory. The results of ^^Fe compound nucleus [5,6] are ideal examples in this situation. 

1 The basic principles to determine p and k from 
intensities of two-step cascades 

Equation (1) contains more unknown parameters than measured the number of inter- 
vals in it. Nevertheless, this type of relation between p and k with J^^ provides limitation 
of the region of their probable values (the simplest example for such type of function is 
giperellipsoid Xl('^2;j)^=const with arbitrary variable number i). The interval of p and 
k possible values is rather narrow only if experimental spectra are decomposed in sole 
primary and sole secondary components according to [3], and ratio r;^^, Vif of primary 
and secondary transitions widths with equal energy and a multipolarity is set in some 
manner. Although the first condition is not satisfied for experimental data from [6], the 
use of algorithm [2] and data from files of ENDS [7] and EGAF [8] permits one to find 
narrow enough intervals of p and k values (with the accuracy {x^ I f < 1) which provide 
reproduction of I^^ in all 27 energy intervals of cascade gamma-transitions. In order to 
get this results, the current k values of iterative process [2] have to be replaced by the 
values which are set by experimental intensities of primary transitions 7.645, 7.279, 6.380, 
6.018, 5.915, 4.217, 3.436 and 3.266 MeV. As usually, in iterative process we used values 
of neutron resonances, low-lying levels densities and total radiative width of s-resonances 
[9]. Level density of ^^Fe at the thermal neutron capture is very small, therefore the 
results of the method used to determine p and k are sensitive enough for parameters and 
decay types of low-lying levels involved in calculation of cascade intensity. Correspond- 
ing data on 10 levels with Ef < 2.113 MeV with spins J'^ = l/2~ — 7/2~ were taken 
from ENDS-file. Radical discrepancy between these results and the p and k values from 
[5,6] has the simplest explanation: reliable data can be extracted from the experiment 
only under condition of execution of basic principles of mathematics and mathematical 
statistics. First of all, should be taken into account specific of errors transfer between the 
measured spectrum and parameters derived from it. From eq. (1), it follows, that in the 
first approach, the error of cascade intensity 
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is the incoherent nonhnear sum of errors 5p and 5k from whole interval of possible energies 
of the cascade intermediate levels. This means, in particular, that the coincidence of the 
experimental and calculated cascade intensities for given p and k in individual intervals of 
nuclear excitation cannot be interpreted (as it was done in [6]) as the proof of their equality 
to unknown values. This fact is illustrated in Fig. 1. Here experimental intensities were 
obtained as superposition of graphic data in figs. 2 and 3 from [6]. They are compared 
with: 

(a) the calculation which uses the results [5] of determination of p and k from the data 
on reaction ^''Fe("^He,"^He,7)^^Fe and 

(b) the typical results of calculation of /^^ for both obtained like in [2] random p and 
k values from the interval of their possible variation, and variant in which level density is 
fixed according to data [5] and only strength functions are varied (see Figs. 2 and 3). 

In case of analysis of the experimental spectrum without its decomposition [3] into 
components corresponding to solely primary and solely secondary transitions, the total 
error of intensity includes in additional error related with registration of cascades with 
secondary transitions from the same energy interval: 51^^^ [Ei) = This leads 

to additional widening of interval of the p and k which reproduce experimental spectra 
due to false solutions. For example, fixation of level density in method [2] according to 
the data [5] results in the sum strength functions which, in the regions ~ 3.2 and ~ 4.4 
MeV, are more than two times larger than the maximum experimental value of k. This 
means that the level density in ^^Fe listed in [5] contains some systematic error. 

1.1 The spectrum of possible functions p, k and their most 
probable values 

Values p and k determined within method [2] contain statistic errors. Apart from 
above classical errors, when calculating the unknown values, the spectra of p and k are 
stretched by the degradation of (1) equation system. This equation system is nonlinear 
and the correlation coefficients absolute values between their quantities are close to unit 
for a given only for some part of p and k values. The remained cases show correlation 
for other pairs of parameters of (1) system. 

It makes impossible situation for degrade linear equations system where interval of 
solution equals infinity. So, in the case considered here, the interval of all possible p and 
k values is always finite. The sign of SI^^{Ei) is different with equal probability at least 
around its zero value) and probability of random localization of iterative process [2] in 
vicinity of point, for example, J2{dl/dp x 6p) = — ^{dl/dk x 5k) is inverse to the value 
Eildl/dp X 5p\ + \dl/dk x 5k\). 

Therefore, one can make conclusion about asymptotically convergence of found solu- 
tions to most probable values, i.e., the average values of ensemble of N realizations of 
iterative process [2] with corresponding standard deviation. Therefore, the p and k func- 
tions estimated in this way differ significantly the from data [5]. Level density (Fig. 2) has 
clearly expressed "step-like" structure as, in practice, in other nuclei studied by means of 
method [2] . The variations of the sum of radiative strength functions of dipole transitions 
(Fig. 3) directly testify to its dependence on the structure of cascade intermediate levels 
in region Bn (if only intensities of cascades [6] with transition energy less than 2.0-2.5 MeV 
do not have large systematic errors). First of all, it should be taken into account that the 
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structure in cascade intensities presented in Fig. 1 in interval 2.5 < Ej < 5.5 MeV cannot 
be reproduced in calculation using level densities and strength functions from [5,6]. In 

correspondence with [3], we could come to a conclusion that data on level densities and 
strength functions listed in [5,6] cannot reproduce total intensity of two-step cascades 
even in principal. In general case, two explanations for the disagreement between results 
of methods [2] and [5,6] could be accepted: 

(a) the presence of serious systematic errors in extraction of p and k from I^^ and 
from the total gamma-spectra corresponding to decay of ^^Fe levels excited in (^He,^He, 7) 
reaction; 

(b) the difference of p and k values reproducing the spectra of these two experiments 
is can result by details of excited levels wave functions structures. Of course, the second 
possibility is to be taken into account only in case of small experimental systematic errors. 

2 Main sources of systematic errors in various meth- 
ods to determine p and k values and any ways of 
their reduction 

So, principal differences in level densities [1,5] and [2] leads to necessity to examine 
data analysis algorithm and to estimate the confident level of results. 

In methods [1,5], common problem at determination of p is the relationship between 
number m = p x AE of levels excited in a given energy interval AE and probability T 
(width r) of measurable reaction product emission. Besides, the corresponding spectra 
can be reproduced by infinite number of functional dependencies p = f{Eex) and T — 
(t){Eex) , the intervals of possible values of level density and emission probability of reaction 
product T can vary from —00 to +00. It means that the type of spectra measured 
with methods like [1,5] always provides the result with lower confident level as compared 
with method [2]. The intensity measured within method [2], to the first approach, is 
inverse to level density and proportional to dk/dE^. Therefore, in all circumstances, this 
provides results with better confidence as compared with methods [1,5]. Again, it should 
be taken into account that, analysis of spectra containing [2] only from events of full energy 
registration of two-step cascades with fixed of their initial and final levels always gives 
considerably better accuracy than analyses the total gamma-ray spectra measured with 
scintillation detector, even though owing to supplementary systematic errors caused by 
subtracting [13] Compton background, and due to higher stability of system with Ge 
detectors. 

In addition, in case [1], determinetion of experimental value p, it required the use 
of theoretical value T of penetration of nucleus surface for evaporated nucleon (or light 
nucleus) . This value should be calculated in frame of nuclear models with error being not 
larger than some tens percents at possible high-frequent sign-changeable variations of T 
relatively to average value. The existence of T gradient directly follows from results of 
modern method [14] to extract p and k from intensities of two-step cascades due to high 
precision of corresponding experiment. As a result, we cannot evaluate confidence of data 
like [1] objectively. Therefore, it is impossible to use these data to confirm [5] results. 

First of all, there is no objective evidence to reject in methods [1,5] the possibihty for 
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compensation of level density increasing (decreasing) in some excitation energy interval 
by the decreasing (increasing) of neutron or gamma-quantum emission probability in 
evaporation or primary gamma-transitions spectrum at the decay or excitation of levels 
close to Eex- This compensation is directly observed [14] for p and k extracted from 
I^^ data. The optical potentials used in method [1], most probably, cannot provide the 
calculation T with required details and errors of p at the level achieved in framework 
of method [14]. The evaluation [15,16] of statistical error and degree of its influence on 
determined in [2] values of p and k showed that method [2] reproduces the peculiarities 
of above parameters for all possible ordinary experimental systematic errors. 

The situation in data analysis from (^He,Q;7) reaction is rather worse. The simple 
analysis [17] of error transfer from primary transition spectra determination to the fit of 
this spectra by "Oslo method" [18] with some p and k values demonstrates that systematic 
errors of the desired parameters increase at least by a factor of several tens at every step. 
At the first step of [17], it is caused by the fact that the intensities of primary transition 
spectrum are much less than intensities of total gamma-spectra (Fig. 4). That is firstly 
related to the low energies of primary gamma-transitions. It is easy to get the lower 
estimation of sum of statistical errors of the primary gamma-transition spectra from the 
narrow energy interval in the vicinity of E^x- This requires one to normalize (per 1 decay) 
all full energy spectra of [17] according to equation: 

J2l'rXE^^E,,. (4) 

The deviation of area of the primary transition spectrum [17] from unity is the part of 
sum of its unknown errors. The maximum of this sum deviation from unit should be very 
small in order to get acceptable systematic errors of p and k. 

At the second step, because of the strong correlation [18] between p and k, the increase 
in its relative errors is practically of the same scale. For example, in order to get error of 
p and k about ~ 50%, it is necessary [19] to determine (and normalize to 1) total gamma- 
spectra at arbitrary Eex with systematic error less than ~ 0.1 — 0.3%. It is possible to 
reject or reduce the systematic errors by using modified variant of method [17] to get 
spectra of primary transitions according to eq. (4) and to determine the best values p 
and k. For example, to use the Gauss method and to correct degrade matrix, this allows 
[20] to realize method of search for the likelihood function maximum within two main 
principles of mathematical statistics: 

(a) to determine supplemental vector of corrections at each iteration by the Jacobean 
matrix without any assumption; 

b) to study explicit form of maximum of p and k likelihood function and to select false 
maxima by simple variation of initial p and k. The method [18] does not fully satisfy two 
above conditions. 

The situation with the confidence level of p and k obtained within method [2] is much 

better. Of course, level density and strength functions derived from experimental I^^ 
data include both ordinary and specific systematic errors. Possible ordinary errors of 
p and k can be easily determined by estimating the propagation error of experimental 
cascade transition intensity to values of p and k. For example, comparison between the 
intensities of the cascade primary transitions from [8] and [21] used for normalization 
/-y^ = F{Ei) and the use of results [22] of extrapolation to zero threshold for distribution 
of random cascade intensities for intermediate levels energies Eex < 0.5 MeV permits 
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precise determination systematic errors for both amplitude and form of F{Ei) function. 
The variation of initial data [2] in limits of their expected errors allows to estimate [20] 

ordinary systematic errors of p and k at presence of nonlinear and ambiguous relation 
between dF{Ei), 6p and Sk. As it was shown in [15], these 6p and dk vahies obtained in this 
manner cannot explain "step-like" structure of p and k [2] by discrepancy between these 
parameters and conventional notions about smooth energy dependence of level density. 

3 The ability to reduce special systematic errors in 
determination oi p k 

The left part of eq. (1) is determined by three unknown functions: the total density 
of the cascade intermediate levels in a given energy interval, sums of radiative strength 
functions of primary and secondary dipole transitions. Strong anti-correlation of these 
parameters causes the phenomenon that total level density with different parity and spin 
(this interval is determined by multipolarity selection rule), and total strength functions 
of El and Ml transitions exactly (x^// << 1) reproducing I^^ vary in narrow enough 
interval. Eventually, this interval for level density with n — +, n — k{El) and k{Ml) 
separately is always wider than for their sums. This statement is true if the primary 
and secondary transitions width ratio is set on the basis of some information. Without 
the above information, the only way to determine p and k in [2] method is the use the 
assumption of equal energy dependencies of k for primary and secondary transitions. 
(Modern nuclear models clearly show incorrectness of this hypothesis at excitation energy 
of about several MeV [23] and, moreover, require [24] detailed accounting for nuclear 
structure at least up to 20 MeV excitation energy). 

The partial compensation of incorrectness of assumption /jP'^«"Y/;;*'ec _ ^onst used in 
[2] is provided by sign-variable deviation of k for secondary transitions at various energy 
with respect to that for primary transitions. It means that the relative variation of 
total radiative width F; of cascade intermediate level in method [2] really is considerably 
less than the relative change the energy dependence of k for secondary transition with 
respect to that for primary transitions. This is observed in all nuclei where experimental 
data permit one to use [14] method to evaluate the general trend of k{E^, E^x) function. 
Furthermore, level density in all nuclei obtained within method [14] has more serious 
discrepancy with the predictions according to [12] than it was earlier established in [2]. 
The effect of change in (^^^""^ / k^'^'^ ^ const) for different energies of decaying level on 
parameters to be determined according to [5] can be stronger because of the used in [18] 
averaging of strength functions for different energies of decaying levels. 

4 Evaluation of confidence level for p and k parame- 
ters and possibility to improve it 

The experimental data on ^^Fe include gamma-ray spectrum following thermal neutron 
radiative capture. These data are suitable, at least, to reveal the difference between the 
experimental and estimated p and k values. Comparison between the experimental and 
calculated spectra for different level densities and radiative strength functions is done in 
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Fig. 5. It confirms conclusion that models like [10-12] cannot be used for precise descrip- 
tion of the compound nucleus cascade gamma-decay. It also testifies to the necessity of 
further analysis of the characters of this process, reduction and more detailed accounting 
for influence of all the sources of systematic errors on the parameters imder search. First 
of all, to reduce error of /^^ = F{Ei) and determine in details energy dependence of func- 
tion k — 4>{E^, Eex) directly from the experiment. This may be done after determination 
of intensities of two-step cascades to final levels with Eex > 0.55„ in new experiment with 
higher quality. Otherwise, one can also to reanalyze experimental data [6]. 

Using numerical method [26] to improve energy resolution together with method [27] 
to determine quantum ordering in resolved in form of peak pairs) cascades and known 
decay scheme [7], one can determine function (1) with small enough systematic error. 
From the population of intermediate levels determined according to [14] from intensities 
of resolved cascades, one can get notion of general tendencies of change in strength func- 
tions as changing energy and involve these data in analysis [2]. This will allow use to 
reject a part of false p and k values obtained in analysis of row cascade intensities and, 
to the first approach, to take into account the change in ratio of intensities of cascade 
secondary transitions to final levels with different excitation energy and deeply different 
wave function structures. 

5 Conclusions 

The analysis (according to the basic principles of mathematics) of the two-step cascade 
intensities following thermal neutron capture in ^^Fe confirms conclusion obtained earlier 
for other nuclei that correspondence between the experimental and calculated [10-12] level 
densities and radiative strength functions within experimental errors cannot be achieved. 
The obtained level density corresponds to general trend observed in [2,14] for other nuclei. 
According to theoretical notions [28], this phenomenon is caused by breaking in this 
nucleus of at least two Cooper pairs of nucleons that results also in change in ratio of 
portions of quasi-particle and vibrational excitations. This phenomenon is not taken 
into account in the models hke [10,12]. Significant increase in strength functions for 
El < 2 MeV completely corresponds [14] to localization of step-like structures in energy 
dependence of level density for secondary transitions of cascades. In totality, the p and 
k values derived from experimental data [6] testify to the effect of nuclear structure on 
these quantities in wide energy ranges and some variation of this effect for different nuclei. 
Perhaps, from position of investigated nucleus with respect to full nucleon shells number. 
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Fig. 1. The two-step cascade intensities in ^^Fe summed over the 0.25 MeV energy 
bins (histogram). Points represent an example of calculation with of p and k obtained 
according to [2]. Calculation with the data from [5,6] is shown by dashed line. The shaded 
area only is used in data analyses [6] . 
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Fig. 2. The total number of intermediate levels of the two-step cascades reproducing 
experimental spectrum shown in Fig. 1. (points with errors). Open points show data [6] 
for both parities and spins 1/2 and 3/2, line represents predictions according to model 
[12]. 
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Fig. 3. The possible strength functions that together with p data from Fig. 2. reproduce 
results of [6] (points with errors). Open points show the data [6]. Curve 1 represents 
k{El) +A;(M1) used to calculate /^^ shown as dashed line in Fig. 1. Line 2 demonstrates 
average k values reproducing together with the fixed p from [5] experimental cascade 
intensities. Dashed line shows predictions of k(El) according to models [10] and [11]. The 
lozenges show the upper estimation of k from the maximum intensities of experimentally 
resolved primary gamma- transitions. 
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Fig. 4. The total gamma-ray spectrum following thermal neutron radiative capture in 
^^Fe (curve 1). Curve 2 represents spectrum of primary transitions for with p and k from 
[5,6] (both spectra are normalized per 1 decay). 
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Fig. 5. The total gamma spectrum summed in the 1 MeV energy intervals (histogram) 
following thermal neutron capture. Points with errors show result of calculation with the 
our data given in figs. 2 and 3. Broken dashed curve shows predictions according to 
models [11,12]. Triangles show calculation with the best [2] values of k and values of p 
from [5]. Open circles are the calculation with p and k from [5,6] only. 
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